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The sequential-reflection model in 
deformed dielectric cavities 

The stationary states of a microlaser are related to 
the decaying quasibound states of the correspond- 
ing passive cavity. These are interpreted classically 
as originating from sequential escape attempts of 
an ensemble of rays obeying a curvature-corrected 
Fresnel formula. Polarization-dependent predic- 
tions of this model, and its limitations for stable 
orbits in partially chaotic systems are discussed. 

As a mechanism for achieving mode confinement, waveguiding by total inter- 
nal reflection is ubiquitous in optics. However, in dielectric microresonators 
where three-dimensionally confined mode volumes are desired, there is al- 
ways leakage because the ray picture, in which Fresnel's formulas describe 
the outcoupling, acquires corrections. Leaky modes corresponding to clas- 
sically confined rays can be found, e.g., in optical fibers as "spiral" modes 
(Poon), or in latterally structured cylindrical VCSELs (Ahn) as well as in 
microdisk lasers (McCall). The classically forbidden loss in such modes is 
analogous to tunneling through an effective potential barrier (Johnson). 

The highest Q is achieved for modes which semiclassically correspond to 
rays almost at grazing incidence. Resonators with a circular cross section 
are a particularly simple realization of this requirement, because they ex- 
hibit whispering- gallery (WG) modes characterized by high intensity in an 
annular region near the surface. However, even Lord Rayleigh who first de- 
scribed the acoustic analog that gave the phenomenon its name, concluded 
(Strutt) that it requires only an everywhere positive curvature, not neces- 
sarily rotational symmetry. A rigorous proof of this is difficult because in 
the short- wavelength limit, this "clinging" of waves to the walls has to carry 
over to the ray picture, in which a generic oval cavity exhibits a transition 
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to chaos (Lazutkin; Nockel, 1996a). Notwithstanding, this problem is funda- 
mental to microresonator design (Nockel, 1997a), because the availability of 
high-Q modes is the foremost selection criterium in an otherwise unbounded 
space of potential resonator shapes (Angelow). Chaos can in fact make WG 
modes more useful, and moreover create other types of modes with desirable 
properties, such as the bowtie pattern whose confocal geometry points the 
way toward the strong-coupling regime in combination with focused emission 
(Gmachl, Morin). 

The robustness of whispering-gallery type intensity patterns in the modes 
of convex resonators extends even to nonlinear media (Harayama). However, 
in that case the distinction to the widely studied phenomenon of vortex 
formation (Weiss) becomes washed out: a WG mode is also a vortex with 
a phase singularity at points of vanishing intensity; for a circular resonator 
where the field is proportional to a Bessel function J m {kr) ~ r m near the 
center r = 0, the vorticity is simply the angular momentum quantum number 
m. Therefore, since our aim is to address the fundamental aspect of the shape 
dependence of high-Q modes in microresonators, we focus here on linear 
media where amplification is taken into account by a negative imaginary 
part of the refractive index n. 

The model considered here can be derived from a homogenous cylinder 
by deforming its cross section and considering only propagation transverse 
to its axis. In this case, TE and TM polarization are decoupled and one has 
to consider only a scalar wave equation 

v 2 i) + fi 2 k 2 i) = o, (i) 

assuming a steady state time dependence so that k is real. Here, n = n — in' 
inside the resonator and n — 1 outside, giving rise to an exterior and interior 
field, ip cxt and ip int , both of which are connected by the proper matching 
conditions at the dielectric interface, depending on polarization. For TM 
modes, ip denotes the electric field, which is parallel to the cylinder axis. In 
this case one finds that ip and its normal derivative are continuous at the 
interface, in analogy to quantum mechanics. 

The system is open because it radiates energy into the environment via 
its modal losses. This openness increases as n —> 1, and the closed-resonator 
limit is approached for n — ► oo. This can be understood from Fresnel's 
formulas which imply total internal reflection for all angles of incidence x 
satisfying sinx > 1/n = sinx c (% c is the critical angle). Equation (1) can be 



recast as 

V 2 ^ + n 2 P^ = 0, (2) 

where n is the real part of n as defined above, and k = k—ik n' / n is a complex 
wavenumber inside the cavity but reduces to k = k outside. If instead of this 
we also had k = k — ikn'/n outside, the solutions of Eq. (2) would be the 
quasibound states of the passive resonator, as they arise when one assumes a 
decaying time dependence oc exp[— ickt — jt], where 7 = ckn'/n. 

For a quasibound (or metastable) state, the field at distances larger than 
~ c/(2-)) from the cavity grows exponentially due to retardation, but within 
this physical range ip ext vanishes as 7 —> 0, so that one can write Vw( r ) ~ 
7C(r). If one expands the dependence of ip- mt and ( on 7 in a Taylor series, 
then to linear order the 7-dependence of (, but not that of ipi nt , can be 
dropped in the full solution. Therefore, the stationary state of the active 
medium and the metastable decaying state are identical to first order in 7 
within an area of order 7~ 2 . 

This approximate equivalence establishes a connection to the study of 
S-matrix poles from which quasibound states arise, see (Schomerus), and to 
dissipation in quantum mechanics (Prigogine, Braun). The recent resurgence 
of interest in these problems is motivated to a significant extent by our lack 
of understanding of the quantum-to-classical transition, in particular in the 
presence of classical chaos. Precisely this constellation is also present in Eq. 
(2) when one considers its short-wavelength limit for the generic case of a 
deformed cavity. 

In the context of laser resonators, there are three main differences to 
previous work on open quantum systems: firstly, we are interested in the 
properties of individual states of an open system, as opposed to a statistical 
ensemble; see also (Casati). Secondly, an important quantity that can be 
studied for such individual states is their emission directionality, which in 
other open systems of chemical or nuclear physics is averaged out. Finally, 
the classical limits of quantum mechanics with smooth potentials and optics 
with discontinuous refractive indices are qualitatively different (Kohler): the 
first yields deterministic Hamiltonian mechanics; the second leads to the 
probabilistic Fresnel formulas which moreover depend on polarization. 

In principle, Eq. (2) can be solved numerically to find the discrete com- 
plex k and the corresponding modes. One approach is based on the Rayleigh 
hypothesis (van den Berg, Barton) which in our implementation for quasi- 
bound states (Nockel, 1996a) assumes that the fields can be expanded in 



cylinder functions as 
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where a polar coordinate system with suitably chosen origin is used. These 
expansions always work inside some circle of convergence for ip int and out- 
side some other circle for ip cx _t, and for a large range of resonator shapes both 
convergence domains contain the dielectric interface where the matching con- 
ditions are imposed to obtain equations for the unknown coefficients A m and 
B m . 

Computational cost can be high here, especially at short wavelengths, and 
hence a semiclassical approximation can lead to simplifications while preserv- 
ing physical insight. The ray picture is a cornerstone of classical optics, but 
its value in the study of open resonators only unfolds when the ray dynam- 
ics is studied in phase space (Nockel, 1994, 1996a, 1996b; Mekis), because 
Fresnel's formulas determine escape probabilities according to the angle of 
incidence \i n °t the position of impact. One can make use of the physical 
information contained in this picture in two ways: Either one starts from 
Eq. (2) and takes a short- wavelength limit (Narimanov); or alternatively, 
one starts from the classical dynamics and makes classical approximations 
that allow one to impose simple quantization conditions and thus make the 
connection to the resonator modes (Nockel, 1997a). The question whether 
these different routes meet "in the middle" is not straightforward because 
the problem of semiclassical quantization in a generic deformed resonator is 
not completely solved as yet, owing to the coexistence of both regular and 
chaotic motion in their classical phase space. 

Among the advantages of the ray-based approach (Nockel, 1997a) are its 
flexibility and computational ease. However, in order for the prescription 
outlined in Ref. (Nockel, 1997a) to correctly describe the limiting case of a 
circular cylinder, one must include the tunneling which in the circle is the only 
loss mechanism. This can be done in the ray picture with a curvature- and 
wavelength dependent "rounding" of Fresnel's formulas which the simulation 
uses at each reflection along a ray path. The idea used in (Nockel, 1997a) 
was to interpret the resonance widths of a circular cylinder in terms of a 
"sequential-tunneling" ansatz: if the intensity of a quasibound state decays 
as exp[— 27 1], this can be interpreted in the ray picture as the result of 



v sequential escape attempts with reflection probability po, where v is the 
number of reflections the ray undergoes during the time t. In a circle of radius 
R, a trajectory characterized by the angle of incidence x nas v = ct/(nL) 
reflections during t (c/n is the speed of light in the passive medium and 
L = 2R cosx is the geometric path length between reflections). Therefore, 
one expects a decay law oc p$ = exp[ct hapo/(nL)]. Comparison with the 
wave result yields 

Po = exp(-2nL'j/c) (4) 

An analytic approximation for 7 in the circle with TM polarization has been 
derived in (Nockel, 1997b), 
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nkR [Jl(kR) + Y*(kR)Y 
Using the semiclassical expression 

m = nkR sin \ (6) 

for the angular momentum (Nockel, 1996a), one then obtains the reflecticity 
in terms of purely classical variables, po(kR, sin x). It reduces to Fresnel's 
formula in the limit of large radius of curvature R, and by construction re- 
produces the width of a mode in the circle if applied locally at each reflection 
in our classical ray model. The latter does not hold uniformly for a similar 
correction derived in (Snyder). 

From the classical limit, it follows that there exists an upper bound on reso- 
nance widths for dielectric cavities with stepped index profiles and "quantum- 
mechanical" continuity conditions on ip, because the reflectivity po, (for po- 
larization perpendicular to the plane of incidence), is bounded away from 
zero. This minimum po,mm will limit the width of resonances in a cavity of 
characteristic size / to ■y ma x = —c lnp 0trnin /(2nl). Smooth index profiles can 
also appear discontinuous on the scale of the wavelength but are eventually 
resolved as k — > 00, allowing arbitrarily small reflectivities at perpendicular 
incidence. 

However, extending these arguments to TE polarization where the electric 
field is in the plane of incidence, we furthermore conclude that a similar up- 
per bound on the widths does not exist even for sharp interfaces. The reason 
is that Fresnel's formula yields zero reflectivity at the Brewster angle \b at 



sinxs = (1 + n 2 )~ 1 / 2 . The normal derivative of ip (which now represents the 
magnetic field) exhibits a jump proportional to n 2 at the dielectric interface 
- reminding us that this is a situation unique to optics. These general con- 
siderations have important implications for microresonator design especially 
at the large n typical for semicondutors, because in that case sinx^ — ► 1/n, 
i.e., the "hole" in the reflectivity for TE polarization approaches Xc for total 
internal reflection from below. Taking tunneling due to finite curvature into 
account as in Eq. (4), the rounded Fresnel formula then exhibits reduced 
reflectivity even for incidence somewhat above Xc- One can approximately 
obtain the TE widths of the circle from 7 = —clm[a]/(nR), 



This is illustrated in Fig. 1, where refractive index and kR are chosen close 
to those of Ref. (Gmachl). The reason is that the quantum- cascade mate- 
rial used there emits preferentially TM polarization, whereas the pioneering 
MQW microdisk lasers with sub-micron thickness permit guiding in the ver- 
tical direction only for TE modes [cf. McCall (1992); there, TE/TM must 
be interchanged to get from the slab- waveguide to our cylinder convention]. 
It is thus important to ascertain whether the identical oval lateral design of 
the quantum- cascade lasers in Ref. (Gmachl) would also permit a microdisk 
laser to operate in TE polarization. 

The lasing mode in Fig. 3D of (Gmachl) was identified as a bowtie-shaped 
pattern corresponding to a periodic ray path with angle of incidence given 
by sinx pa 1/n, i.e., directly at the critical angle. That this mode provides 
high Q can be seen by comparing to 7 in Fig. 1: assuming that the width 7^ 
of a bowtie mode results from the sequential application of po as determined 
for the circle, the argument leading to Eq. (4) implies that 7^ ~ 7 L/l where 
L/l pa 1.13 is the ratio of the classical path lengths between reflections in 
the WG orbit and bowtie, respectively. One sees that the TM line intersects 
the critical angle (corresponding to m — kR) at a much smaller width than 
the TE curve, and this Q-spoiling due to Brewster transmission is borne 
out by the actual TE resonances as well. This leads to the prediction that 
conventional microdisk lasers with a shape designed to yield a bowtie pattern 
just at Xc as in Fig. 3D of (Gmachl) will not lase. 

These ray arguments are known to yield large deviations from the true res- 
onance widths when the modes under consideration are quantized on stable 
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Figure 1: Exact resonance widths of a circle (radius R, refractive index 
n = 3.29) at 35 < kR < 35.5 versus angular momentum m. Dashed line 
indicates \ c in Fresnel's law, using Eq. (6). Also shown are the TM widths. 
The Brewster angle is at m ~ 34 and causes a peak in the TE widths. 

phase-space domains in a partially chaotic system, cf. Ref. (Nockel, 1997a) 
where this was attributed to chaos-assisted tunneling. The latter yields en- 
hanced outcoupling and hence the true widths are underestimated by the 
sequential ray picture. Therefore, the above Q-spoiling for TE modes is 
not counteracted by a correction of this nature. The prediction of an upper 
bound for TM widths is also not affected by chaos-assisted tunneling because 
it cannot be faster than the fastest classical process, which in turn is limited 
by po at sinx = 0. Beyond this, however, quantitative widths for stable-orbit 
modes in mixed phase spaces are not provided by the ray model. 

The disagreement is illustrated in Fig. 2 for a bowtie mode similar to the 
ones studied in (Narimanov), as a function of n, but at a deformation of e = 
0.16 [defined as in (Gmachl)] and nkR m 119.8. Since nk is the wavenumber 
inside the resonator, it should remain approximately independent of n as 
long as the outcoupling can be taken into account in the form of a boundary 
phase shift intermediate between Dirichlet and Neumann. Indeed, for the 
state shown in Fig. 2, the change in nkR in the plotted range of n is only 
« 0.2. The length scale R here is the radius of curvature at the points of 
reflection. At small n where Xc is larger than the angle of incidence of the 
bowtie, escape is classically allowed in Fresnel's formula and hence curvature 
corrections are unimportant. At n > 3, the tunneling correction in Eq. (5) 
does improve on the classical Fresnel prediction (7 = 0) but clearly still 
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Figure 2: Width of a TM bowtie mode vs. refractive index, from numerical 
and ray calculations. 

underestimates the true width. As tunneling in general is definable only 
with respect to a classical expectation, we could again label the discrepancy 
as chaos-assisted tunneling. However, a semiclassical theory starting from 
Eq. (2) which reproduces the exact behavior in Fig. 2 very well (Narimanov) 
can shed more light on the physics of the phenomenon. 
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